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Abstract 

The dual of an arbitrary D-dimensional nonabelian lattice gauge theory, obtained 
after character expansion and integration over the gauge group, is shown to be a 
local lattice theory in the eigenspace of the Casimir operators. For D < 4 we also 
provide the explicit form of the action as a product of character expansion coefficients 
and Racah coefficients. The representation can be used to facilitate strong coupling 
expansions. Furthermore, the possibility of simulations, at weak coupling, in the dual 
representation, is also discussed. 
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Ever since their invention by Wilson |jT[, lattice gauge theories have inspired research in 
a multitude of directions. One of the most interesting of these has been the investigation of 
duality properties. The study of duality properties of various spin models has a long history, 
but more recently duals of compact U{1) lattice gauge theories were also found 

The first step of finding duals of lattice gauge theories is always the character expansion 
of the Wilson-action. The second step is integration over the gauge degrees of freedom. In 
the case of U{1) theories these steps lead to integer valued fields satisfying constraints. The 
constraints have the form of linear equations for the eigenvalues of Casimir operators. For 
f/(l) theories these are additive quantum numbers. 

Suppose now that the dual picture is to be used for simulations. Then either the additive 
quantum numbers involved in the constraints have to be summed over, or the constraints 
(Bianchi identities) have to be resolved in terms of fields with no constraints. The first of 
these routes is not feasible for U{1) gauge theories: it would be tantamount to the exact 
solution of the theory. However, the second route has been successful 0. 

When a similar procedure is applied to nonabelian lattice gauge theories the above two 
steps also lead to local discrete valued models. The partition function of such a model is 
formed from a sum over products of the character expansion coefficients, depending on the 
eigenvalue of the quadratic Casimir operator and of Clebsch-Gordan coefficients depending on 
two kinds of discrete valued fields. Fields belonging to the first group take values that are the 
eigenvalues of diagonal elements of the Lie algebra, labeling states in a given representation. 
We will denote the collection of them for a particular group element by m, n, M, N, ^, or 
u. These are additive quantum numbers and they all satisfy constraints in the form of 
linear equations. Either these constraints have to be resolved, or the summation over the 
corresponding fields has to be performed. The second group of discrete valued fields take 
values that are the eigenvalues of Casimir operators. They will be labeled by letters j, J, 
and /. They satisfy triangular and tetrahedral inequalities. As will be pointed out later, 
such constraints would allow simulations respecting detailed balance. For short we will refer 
to these fields as magnetic quantum number and angular momentum, respectively, though 
these names are truly appropriate only if the gauge group is SU{2). 

The purpose of this letter is to investigate the possibility of performing summations over 
magnetic quantum numbers. Of course, this can be done in a formal manner, but the crucial 
question is whether the resulting effective theory of angular momentum valued fields is local 
or not. We will answer this question in the affirmative for arbitrary compact gauge group 
and arbitrary space-time dimension D. In fact, the statement is fairly obvious for the more or 
less trivial case D = 2, and it has recently been proved for D = 3 hj several authors 0. 

A lattice gauge theory in the fundamental representation (a generalization to other rep- 
resentations is possible) has the form 



where Xf denotes the character of the group in the fundamental representation depending 
on the plaquette element, to be described later, and / dAi{r) represents integration over the 



Z = n II dA{r) n exp{/?x/[(f/.fc(r)]} 
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Haar measure. The product index i represents the axes of the lattice, while the (ik) pair, in 
conjunction with r, represents a particular plaquette (2 dimensional simplex) of the lattice. 
The exponentials can be expanded in a series of characters to give 

exp{/?x/[(t/«fc(r)]} = EdjXAUrk{r)]c,{p) (2) 

j 

where dj is the dimension of the representation labeled by j. 

Now the characters are traces of the representation matrix of the plaquette element. The 
plaquette element can be written as the product of the link elements at the boundary of the 
plaquette. The link element Ui{r) is associated with a link originating at r and extending one 
lattice unit into the direction of the 2th axis. The element Ui{r) is also associated with the 
positive direction, while its adjoint with the negative direction. Setting r = 0, for simplicity, 
we have 

xAUikiO)] = Trj[uH0)Ukmi(ek)Ulm, (3) 

where is a unit lattice vector in the positive direction along the kth axis. 

The dual of the lattice consists of D dimensional hypercubes centered at the vertices of 
the original lattice. The boundary of a hypercube consists of 2D codimension 1 cubes. The 
dual of a link of the original lattice is such a cube, on which the link variables [/,; live. The 
duals of plaquettes are codimension 2 simplexes (for D=4, plaquettes themselves). They 
are the boundaries of the D — 1 dimensional cubes that themselves form the boundary of 
the hypercube. Duals of plaquettes attached to a single link of the original lattice form the 
bondary of such a, D — 1 dimensional simplex. There are 2D{D — 1) such simplexes attached 
to each hypercube. 

Our discussion starts with a simple observation: 

1 ) Every pair of adjacent link variables of a plaquette define a unique hypercube. 
The four possible pairs one can form on a loop define four adjacent hypercubes. 

The proof is simple. Two adjacent links of the original lattice intersect in a lattice point. 
Let us associate this point with the pair of adjacent link variables. A point of the original 
lattice is, however, dual to a hypercube (centered at the point). The corners of a plaquette 
define four separate, but adjacent hypercubes. Their mutual boundaries are the duals of the 
link, codimension 1 simplexes. 

Now saturate each product of link matrices on the plaquette by a complete system of 
states. Each complete system can be labeled by (j, m), where m represents the collection of 
quantum numbers labeling states in representation j. Then it follows that 

2) Every magnetic quantum number, m is associated with a unique hypercube. 

This is true because every m is associated with a pair of link matrices, but according to 1) 
every pair of adjacent link matrices is associated with a hypercube. 
The plaquette contribution can then be written as 

X,[[/..(0)] = i^^,„,(t//)i?^,^3(f/.)i^^3m43(f/^)^^,„.,(C^fe^), (4) 
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where the arguments of the hnk operators were dropped for simphcity. They are the same 
as in d^). Each of the four indices is associated with an adjacent hypercube, according 
to 1) and 2). 

Fig. 1. shows a plaquette in the ik plane. It also shows the magnetic quantum num- 
bers labeling systems saturating adjacent products of ^/-operators. The magnetic quantum 
numbers are associated with dual hypercubes centered at corners of the plaquette. 

It is worth making a comment at this point. Lemma 1 fails for globally symmetric 
nonabelian spin models. There is no natural way of associating the two hypercubes, con- 
nected by a link, with one or the other of the set of variables nii and m2 of the product 
-^mim2(^^)-^m2mi(^')- Indeed, their duals are not local in the space of Casimir operators. 
The results of this paper would fail for nonabelian spin models. 

Furthermore, it is easy to see that 

3) Every link matrix Ui is associated with the same hypercube on its left (and also 
on its right) in every plaquette it enters. 

The proof follows from the definition of link variables. Naturally, the association is 
reversed for the matrix Uj . Then we obtain 

4) For every link variable, U, the left subscript (right subscript) of each rotation 
function Dl^,^{U) it enters is associated with the same hypercube. 

4) is a simple consequence of 2) and 3). 

Let us prove now the most important result of this letter: 

5) Every Clebsch-Gordan coefficient, obtained after combining rotation functions 
and integrating over the gauge variables, depends on magnetic quantum numbers 
associated with a single dual hypercube only. The angular momentum variables 
of the Clebsch-Gordan coefficients are shared by neighboring hypercubes only. 

The proof starts with the use of the addition theorem for representation functions to 
reduce the number of representation functions {2D — 2) depending on the particular link U. 



J 



(5) 

where the Wigner's three-j symbols have been used rather then Clebsch-Gordan coefficients. 

When a rotation function with argument U is combined with another one with argument 
W , then the relation 

DUU^) = (6) 

should be first used. 

Two important observations concerning should be made at this point. First, the 
dependence on indices m and n, associated with different hypercubes, factorizes in the two 
three-j symbols. Second, the new rotation function appearing in the addition theorem also 
satisfies 4), because M = nii + m2 and N = ni + n2. 
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Subsequent applications of the addition theorem reduce the product of 2 (D — 1) rotation 
functions to a single rotation function D-^^{U) and a product of 2{2D — 3) three-j symbols, 
2D — 3 of which depend on 2(1^ — 1) the quantum numbers only, which are associated 
with one of the hypercube and the same number of them depend on quantum numbers 
only, which are associated with the other hypercube. Finally, integration over the group 
implies that J = /i = z/ = Oin the last rotation function. This, in turn, implies that the 
last two three-j symbols that contain the quantum numbers (J, /x) and (J, z^), respectively, 
turn into Kronecker deltas for two pairs of angular momentum vectors. That leaves us with 
the product of 2{D — 2) three-j symbols dependent on quantum numbers rrii only and the 
same number of symbols dependent on the quantum numbers rii only. The dependence on 
additive quantum numbers labeling states is completely factorized. 

Notice now that a similar construction can be performed on every one of the 2D links 
emanating from a given lattice point. Altogether, we will have AD{D — 2) three-j symbols 
dependent on rrij quantum numbers associated with a given hypercube only. Furthermore, 
the angular momenta in these three-j symbols are shared by the two hypercubes joining in 
the D — 1 dimensional simplex, dual to the appropriate link variable, only. This completes 
the proof of the theorem. 

Now the final form of the main result of the paper can be spelled out as 

6) Summation over magnetic quantum numbers results in a discrete local field 
theory in the eigenvalues of Casimir operators. 

Summation over the magnetic quantum numbers results in factors associated with hy- 
percubes. The factor depends on angular momenta corresponding to the 2D codimension 1 
simplexes of the hypercube and additional 2D{D — 2) angular momenta obtained at repeated 
applications of addition theorem (|^), also associated with the same hypercube. Thus, the 
range of the interactions is the size of a single hypercube. 

Let us investigate now the simplest cases D = 2,3, and 4. Although D = 2 lattice gauge 
theories are trivial, still the application of our results to them is instructive. Also the algebra 
is simple enough so that it can be written out in detail. There are D = 2 link variables Ui 
and U2 running out of an arbitrary point of the original lattice and D = 2 variables Ui and 
U2 running into it. The 2D{D — 1) = 4 plaquettes containing these variables are 

TriU^Ui...), Tt{uIU2...), Tt{uIuI..), Tr(/7it/j...), (7) 

where link operators connecting other dual plaquettes (hypercubes) have been omitted from 
the traces. 

After inserting complete systems of states one obtains the following product of rotation 
functions: 

J:m,,m,,m,,m,, D^.^LjU^) D^^l^^JU^) D^^^^JU^ 

where indices associated by other dual plaquettes have been omitted. 
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Integration over the rotation functions is simple, because there are only two functions 
for each rotation, each giving a Kronecker delta for the angular momenta and magnetic 
quantum numbers as well. Thus, the summation over magnetic quantum numbers results 
in a factor 2j + 1 only, where j = ju = ■■■ = ji2- The factor corresponding to a dual 
plaquette is proportional to a Kronecker deltas for all the angular momenta involved and a 
trivial j-dependent multiplier. Since the angular momenta are shared between neighboring 
dual plaquettes, the partition function becomes diagonal in angular momentum. Thus, the 
partition function is Z]j(2j + l)[cj(/5)]^, where V is the volume of the lattice. 

The D = 3 case |^ is more complicated. There are 2D = 6 link (dual plaquette) 
variables, and 2D{D — 1) = 12 plaquettes (dual links) involved. Each link variable appears 
in four different plaquettes. Each appearance contributes to the expression by a rotation 
function, as in (|[). The four rotation functions can be combined pairwise using (^ and (D 
giving two three-j symbols for the dual plaquette in question and two rotation functions. 
Integration over the group space results in the identification of the new angular momenta and 
magnetic quantum numbers in the three-j symbols. Before we write down the appropriate 
expressions obtained after this procedure we introduce a concise notation, due to Wigner 0. 
Magnetic quantum numbers will be omitted from three-j symbols, with the understanding 
that repeated angular momenta in products imply summation over the corresponding mag- 
netic quantum numbers. Since we restrict our discussion to three-j symbols associated with 
a single hypercube, every angular momentum will be uniquely associated with a magnetic 
quantum number.^ 

Labeling angular momenta as in (H) we obtain the following product of three-j symbols: 

(jl3il2-^l) (il2jl3-^l) (il2j23>^2) (ji2i23-^2) UlshaJs) {hsj 13 -h) 
X (Ji2il3^i) (il3Jl2^l) (il2 ^23^2) (j23il2^2) (ilsJSS-^s) {kdl^J^), (9) 

This is the complete expression for the contribution of a dual cube, apart from the multiplier 
n(2Ji + 1), shared by neighboring dual cubes. 
Using the relations 

(JU2J3)(J2J3^1)(J3J1J2) = ( f ^']{JM) (10) 

Jl J2 33 



and 



UlJ2J3){j2j3Jl)U3jlJ2)iJlJ2J3) = \f f M (H) 

Jl J2 J3 



reduces to 



J-l J 2 J3 1 I Jl J-2 ^3 11 Jl J2 J-3\\ J-l J-2 J-3]\ J-l J-2 J-3 
j23 jl3 jl2 J 1 i23 jl3 jl2 J 1 ^23 jl3 jl2 ] \ Jl J2 J3 ] \ M jl3 jl2 



(12) 



^ As Wigner points it out |g] contravariant and covariant components should be used in the three-j 
symbols, which differ in sign. Contravariant component are always contracted with covariant ones. Our final 
result is not affected by this complication, since we intend to sum over all magnetic numbers. Thus, for the 
purpose of simplifying notations we will not distinguish contravariant and covariant components. 
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(p!^) has a simple geometrical interpretation if the gauge group is SU (2) [|| [Q. Since angular 
momenta are attached to plaquettes on the original lattice, and links (vectors) on the dual 
lattice, each six-j symbol represents a tetrahedron. The five tetrahedrons corresponding to 
the five six-j symbols form a cube-like object with triangular faces. The angular momenta 
Ji and J'i correspond to face diagonals of the dual cube. All angular momenta are shared by 
neighboring dual cubes. The angular momenta jik are shared by three other, neighboring 
cubes: the ones displaced by ej,efc, and by ej + e^. Barred indices correspond to negative 
directions. Each of the angular momenta Ji is shared only by one other cube, the one 
displaced by vector e^. The collection of tetrahedra cannot always be embedded in fiat three 
dimensional space 0. In fact, if the coefficients Cj{(3) of (H) are omitted then the partition 
function reduces to the discrete model of Regge for three dimensional gravity. |^ |^ 

Let us turn now to the physically most interesting D = A case. Integration over the gauge 
degrees of freedom leads to the following invariant sum over products of three-j symbols: 

(jl2i244^) (jl2il4^D 0^4^244^) (^34^23^1'^) (j24i23 4"^) (j24i34^f) 

X {jiiiuJl'^) (ji2ii4<^P) (ji4i24'4^) (j'siisi'^^) (issisi'^l'^) {h^hiJl^) 

/ . . -rT2\ / . . -r'i.As I ■ ■ -A'X\ I . ■ t14\ / . . T^ZlN / ■ • tI^N 



X 



X 



X 



X 



(Ji3j344') (Ji3jl4 ) (Ji4j344') (il3j34 4') (jl3Jl4 ) (jl4j34 ^1') 
(Jl2i23 J-2 ) (Jl2il3 ) (Jl3i23 jf) {3X2321 J^2 ) {3X2311 jf ) (ji3i23 ) 



\ L 1 1 f \ L Z Z / V O Hi O J \ ^ 4 ^ 

I 724 t34 t23n / t14 t34 t_13w t24 t_14 t12w t_12 t13 t2 
Wl "^X "^1 JW2 "^2 "^2 JW3 "^3 "^3 ^"^4 "^4 "^4 



12 7-13 T23^ 
4 ^4 ; 

13 t23n 



(13) 



Again the symbols j correspond to angular momenta carried by plaquettes, while symbols J 
represent angular momenta obtained when the j-s are added. Both of these angular momenta 
are shared by neighboring hypercubes, exactly the same manner as for D = 3. Note that 
the subscript of symbols J indicates the direction of the neighboring hypercube sharing the 
corresponding angular momentum. 
Using relation (|10|), (|TT]), and g] 



{3\323){hH3) 



.1)2.3 ^(2/ + 1 




{kH3'){3z323') 



(14) 



one can express product (|T3|) by means of invariant, six-j symbols. We obtain 
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^ { Tl2 T 723 } { Tl4 Tl3 7 "i t24 723 734 (15) 





J5 

where the subscripts of the Jf'' angular momenta, appearing in (|l^) have been omitted. The 
angular momenta Ji,...,Ji appear in the process of reduction of (P^. They are not shared 
by neighboring hypercubes. 

Finally, we examine the possibility of simulations, using the dual representation. In 
principle, such a possibility is exciting, because simulating integer valued theories without 
the need of tedious matrix multiplications makes computations faster. There are several 
issues, however, which have to be resolved before simulations can be attempted. 

The first issue is detailed balance. Suppose an angular momentum j is updated. Suppose 
A is the range of the allowed random change, i.e. the randomly proposed new value of j 
satisfies j — A < j' < j + A that is independent of the inequalities that j satisfies. Then there 
are two possibilities. Either / also satisfies all the required inequahties, or it does not. In 
the former case the probability of the j j' transition is the same as that of the transition 
j' — > J. In the latter case the two probabilities are still equal, namely they are zero. Thus 
for such an update procedure the detailed balance condition is satisfied. 

The second issue is the complicated form of six-j symbols. Note, however, that the aim 
is to perform simulations in the weak coupling regime, relevant in the continuum limit. At 
weak coupling the average angular momentum is large. In fact, j ~ ^/^. For large values of 
the angular momentum the semiclassical limit of six-j symbols can be used. Wigner has 
shown that in that limit six-j symbols are given in an average sense by 

/ ^2 ^3 Y 1 .^g^ 

1 Jl J2 Js J "47r|(jiXj2)-J3|' ^ ' 

i.e by the inverse of 247r times the volume, V, of the tetrahedron. The volume can be 
expressed by the edges of the tetrahedron using Cayley's formula 
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(17) 



The volume has square root type zeros at the edge of the allowed region, thus, the six-j 
symbols have mild, integrable singularities as the tetrahedron becomes degenerate. 

One may also consider using continuous angular momentum variables. This would be 
approximately equivalent to using a noncompact gauge group. It should not alter results in 
the weak coupling limit. 

One very important issue we have not yet resolved is the positivity of the integrand of 
the functional integral in j-representation. Indeed, (p!6| ) gives an expression for the square 
of six-j symbols only. The six-j symbols themselves oscillate. The phase of oscillation was 
found by Ponzano and Regge |0| In the semiclassical limit, summations over angular 
momenta are dominated by values which make all phases of oscillations stationary. The form 
of the resulting field theory is complicated (at least in the D = A case) and will be dealt 
with in a future publication. 
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Figure Captions 

Fig.l. A plaquette in the ik plane, mj represent magnetic quantum numbers associated with 
dual hypercubes. 
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